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Abstract 

We study inflationary solutions in the M-theory. Including the fourth-order curvature 
correction terms, we find three generalized de Sitter solutions, in which our 3-space ex- 
pands exponentially. Taking one of the solutions, we propose an inflationary scenario of 
the early universe. This provides us a natural explanation for large extra dimensions in a 
brane world, and suggests some connection between the 60 e-folding expansion of inflation 
and TeV gravity based on the large extra dimensions. 
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1 Introduction 



Recently there has been renewed interest in obtaining accelerating universe from higher- 
dimensional gravitational theories. The study of this subject gets its urgency from the 
recent discovery of the accelerated expansion of the present universe as well as the con- 
firmation of the existence of the early infiationary epoch [1]. Though it is not difficult to 
construct cosmological models with these features, it is desirable to derive such a model 
from the fundamental theories of particle physics that incorporate gravity. The most 
promising candidate for such theories are the ten-dimensional superstrings or eleven- 
dimensional M-theory, which are hoped to give models of accelerated expansion of the 
universe upon compactification to four dimensions. 

It has been shown that models with certain period of accelerated expansion can be 
obtained from the higher-dimensional vacuum Einstein equation if one takes the internal 
space hyperbolic and its size depending on time [2]. It has been shown [3] that this class 
of models is obtained from what are known as S-branes [4, 5] in the limit of vanishing 
flux of three- form fields (see also [6] ) . It is found that this wider class of solutions gives 
accelerating universes for internal spaces with zero and positive curvatures as well if the 
fiux is introduced. Further discussion of this class of models has been given in refs. [7]- [9]. 

It turns out, however, that the model thus obtained does not give enough e-folds 
necessary to explain the cosmological problems such as horizon and fiatness problems [3, 7]. 
Some efforts to overcome this problem were made in the present framework in ref. [9]. 

We note that the scale when the acceleration occurs in this type of models is basically 
governed by the Planck scale in the higher (ten or eleven) dimensions. With phenomena 
at such high energy, it is expected that we cannot ignore higher order corrections in 
the theories at least in the early universe. In fact there are terms of higher orders in 
the curvature to the lowest effective supergravity action coming from superstrings or M- 
theory. In four dimensions, many studies have been done with such correction terms [10] . 

The cosmological models in higher dimensions were also studied intensively in the 
80 's by many authors [11, 12]. Among them, the model with Gauss-Bonnet term is very 
interesting for the above reason [12]. It was shown that there are two exponentially 
expanding solutions, which may be called generalized de Sitter solutions since the size 
of the internal space depends on time (otherwise there is no solution of this type). In 
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both solutions, the external space inflates, while the internal space shrinks exponentially. 
(There are also two time-reversed solutions, i.e. the external space shrinks exponentially 
but the internal space inflates.) One solution is stable and the other is unstable. Since 
the present universe is not in the phase of de Sitter expansion with this energy scale, we 
cannot use the stable solution for a realistic universe. If we adopt the unstable solution, 
on the other hand, we may not find sufficient inflation unless we fine-tune the initial 
values. The higher-order curvature terms called Lovelock gravity were also considered in 
higher-dimensional cosmology [13]. 

However, most of the work considered powers of scalar curvature or Lovelock gravity, 
which are not the types of corrections arising in type II superstrings or M-theory. In 
particular, it is known that the coefficient of the Gauss-Bonnet terms vanishes and the 
first higher order corrections start with terms (one is the fourth order Lovelock gravity 
and the other contains higher derivatives) [14]. The purpose of this paper is to examine 
how these corrections in the fundamental theories modify the above cosmological models 
and whether we can get interesting cosmological scenario with large e-folds. We focus on 
the solutions to the vacuum Einstein equations with these higher order corrections in this 
paper since the basic feature can be obtained in this simple setting. 

We find interesting models with power law as well as exponential expansions. These 
seem to give enough e-folds, and the solutions suggest that the fundamental Planck scale 
is rather small as < 10^ TeV and the size of internal space grows rather large with the 
scale < 10 TeV"^. These may provide the models of large extra dimensions discussed 
by Arkani-Hamed et al. [15]. We have found similar results for several superstrings and 
M-theory, but here we present only the results of exponential expansions for M-theory, 
leaving other details to a separate paper [16]. 

2 Vacuum Einstein Equations with corrections 

We consider the low-energy effective action for the M-theory: 



S = 



-S'eh + 'S'4, 




where 




(2) 
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= ?V / d}^x^Ujo + 2l3E^ . (3) 



P„ al3^HlVi...^l4,V4, r>pl(Tl . . . PP40-4 (A \ 

8 — 2^ X 3! ta^^pjCTl...p40-4-"' /ili^l ^41^4 ) l^/f 

Jo — R^^'^'^Ra/xi^pRx^'^'^R^ paK + -^R^'^^'' RafBpvRx'^'^" R'^ pan- (5) 

Here we have dropped contributions from forms, nf^ is an eleven-dimensional (IID) grav- 
itational constant, and we leave the coefficients /3 and 7 free although we know them for 
the M-theory [14] as 

^ 32 X X (27r)4' ' 3 X 24 X (27r)4 ' ^ ' 

where T2 — {2t^'^ / k\^^/^ is the membrane tension. Type II superstring has the same 
couplings in 10 dimensions, so we can discuss this case if we keep the dilaton field constant, 
but we consider IID theory in this paper. Here we should note that contributions of the 
Ricci tensor R^i, and scalar curvature R are not included in the fourth-order corrections (3) 
because these terms are not uniquely fixed. 

Since we are interested in a cosmological time-dependent solution, we take the metric 
of our spacetime as 

i=l a=5 

where we assume that the external 3-spacc and the internal 7-spacc arc flat. Taking 
variation of the action with respect to Mq, Ui, and wc obtain three basic equations, 
whose explicit forms will be given in a forthcoming paper [16]. 

2.1 Generalized de Sitter Solutions 

In cosmology, de Sitter inflationary expanding spacetime is the most important solution 
in the early universe. Hence, let us flrst look for such solutions. Assuming the metric 
form of a generalized de Sitter spacetime as 

Wo = 0, Ui = nt, U2 = I't, (8) 
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where /j, and v are some constants, we obtain three algebraic equations: 



+ Tjiu + z/^ 



-77 



+47(3// + 7i/) 6//^ + 42i/^ + 7/iV^(^ + + i/^ + = 0, 



(9) 



3/^2 + + 28z/2 + 20160/3z/^ 
+37 



27/x(3/i + 7z/) 48/i^ + 7z/2(3/i^ + 2/iz/ + z/^)(/i2 + z/^ + /xz/) 

+ 14//V2(2// + i/)(3// + i/) + 42i/^(// + 2i/)(// + u) 



+27/x(3/x + 7z/)^ 12/i^ + 7z/2(/i + z/)(/i2 + z/^ + /xz/) 



0, 



(10) 



2/x^ + 611U + 7i/^ + 2880/?//i/^ 
+7 



15//^ + 46/xV + 38//i/^ + 6i/'' 
12/ + 7/i/^(/ + i/^ + //i/)^ + 168u^ + 7//V^(2// + + 2l/i/^(// + 21/)^] 
27z/(3/i + 7z/) 96z/^ + + 2/xz/ + 3z/2) (^^2 + ^,2 ^ ^j,) 

+ 2/(2// + i/)(// + I/) + 6/i/2(// + 2i/)(// + 3z/)] 



+27z/(3Ai + 7z/)^ 12z/^ + /(/x + z/)(/ + z/^ + /xz/) 



0. 



Since these equations are very comphcated, we have solved them numerically. If 7 
does not vanish, rescaling /3, 7, and u as 



P = /5/I7I , 7 = 7/I7I (= 1 or - 1) , /t = Ai|7|^/^ , and z> = u\-f\ 



1/6 



(12) 



we can always set 7 to —1 if it is negative (or 1 if positive). We also have to rescale time 
coordinate as i — \^\~^/H. The typical dynamical time scale is then given by \^\^^^ ~ 
0.181818mjf/, where mn = ^n^^ is the fundamental Planck scale. After this scaling, 
we have only one free parameter /3. In Fig. 1, we depict numerical solutions Nj(/ij, z/j) 
(i = 1 ~ 5) with p,i > with respect to P for the case of 7 < 0. We note that there 
are always time-reversed solutions N^(/, j^l) (i = 1 ~ 5) obtained by (/i^, i>l) — —{p,i, Vi) 
which are not shown explicitly. We find that M-theory (7 = —1, /3 = /^^ = — 1/(3 x 2^) fa 
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Figure 1: Five generalized de Sitter solutions Nj(/ij,z/j) (i = 1 ~ 5) with /ij > with 
respect to /? for 7 < 0. Each pair of points (/ij, z/j) with the same value of /? gives one 
solution. Number of solutions changes with the value of (3. There is another set of 
time-reversed solutions N- with /ij < 0. 

—0.0052083) has three solutions 

N2(/i2,i>2) = (0.45413,0.45413), 

N3(/i3,i>3) = (0.79802,0.10781), (13) 

N4(/24,i>4) = (0.50754,0.43025). 

If 7 > 0, we find three solutions Pj(i = 1 ~ 3) for /9 < 0, while just one P4 for /3 > 0. 
We depict these solutions in Fig. 2 (for 7 = 1). There is no solution for /? = 0. In Table 
1, we summarize their properties. 

If 7 = and /3 7^ 0, we can always set /3 to —1 if it is negative (or 1 if positive), by 
rescaling /?, /i and z/ as /? = (= 1 or — 1) ,/2 = , and v = . We then 

find two solutions [{fi, v) = (0.22046, -0.28771), (0.19168, 0.19168)] for (3 = -1, while just 
one [(/X, z/) = (0.04467, -0.45111)] for (3 = 1. There is no solution for /5 = 7 = 0. 

2.2 Stability 

Since the solutions obtained above correspond to fixed points in our dynamical system, 
we have to analyze their stabilities in order to see which solutions are more generic. 
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Figure 2: Four generalized de Sitter solutions Pj(/ij,i>j) (i = 1 ~ 4) with jli > with 
respect to /5 for 7 > 0. Each pair of points (/ij, z/j) with the same value of (3 gives one 
solution. Number of solutions changes with the value of (3. There is another set of 
time-reversed solutions with /i^ < 0. No solution exists for 7 > and (3 = 0. 

We have performed a linear perturbation analysis around those fixed points. Setting 
dui/dt = fxi + AiC"^ and du2/dt = z/j + BiC"^ , where l^jl, -C 1, we write down the 
perturbation equations. There are five modes (cr = aa\ a = 1, 2, • • • , 5) because the basic 
equations for iii and U2 are two third-order differential equations plus one constraint which 
is second order. We show the results in Tables 1 and 2. In Table 1, we just give the number 
of stable and unstable modes. For example, there are one stable and four unstable modes 
for the solution Ni(/2i, i>i). Hence this solution may not be generic because we have many 
unstable modes. The M-theory has three solutions (13). Two solutions (N2 and N3) have 
four stable and one unstable modes. The third solution (N4) has five stable modes, which 
means that this solution is stable against hnear perturbations (see Table 2). 

From this stabihty analysis, one may conclude that the solution N4 is most preferable 
spacetime in this model. We have inflationary expansion not only in 3D external space 
but also in 7D internal space. The expansion rates in both spaces (/i, u) are, however, 
almost the same. Hence, we would find our present world in which scales of two spaces are 
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Table 1: Generalized de Sitter solutions Nj(/ij,i>j) (i = 1 ~ 5) and Pi{p,i,i>i) (i = 1 ~ 4) 
with jli > for various values of $ in the cases of 7 = —1 and of 7 = 1, respectively. 
The solution N4 is the partner of the solution N4 with p,4 < 0, which is found when it 
is extended to the range of ;9 > 0.025668. Five eigen modes for linear perturbations are 
also shown. (ms,nu) means that there are m stable modes and n unstable modes. The 
solution has many stable modes if its 10-volume expansion rate (3/i + 7u) is positive. 
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solution 


property 


range 


stability 


3/i + 7i} 






/7i > D > 771 


^ < _0 04342 


fls 4u) 






N2 


p,2^ i>2> 


(3 < -0.00418 


(4s, lu) 


+ 






-0.00418 <P < 0.0004464 


(5s,0u) 


+ 


7 = -l 


N3 


fi3> i>3> 


-0.202867 < P< 2.141 


(4s,lu) 


+ 




N4 


/i4 > i>4 > 


-0.202867 <P< -0.00418 


(5s,0u) 


+ 






Z>4 > /i4 > 


-0.00418 <P< 0.025668 


(4s, lu) 


+ 






/i^ > > i>4 


P > 0.025668 


(ls,4u) 






N5 


/is > > i>5 


0.025668 < /3 < 0.09247 


(ls,4u) 






Pi 


jli > > ui 




(ls,4u) 




7=1 


P2 


p,2 ^ i>2> 


/3 < 


(4s,lu) 


+ 


P3 


/ii > > z>i 




(4s,lu) 


+ 




P4 


/ii > > z/i 


/3 > 


(ls,4u) 





not so different. This solution also predicts that inflation never ends because it is stable. 
For a realistic cosmological model, the solution must be unstable because inflation should 
end. On the other hand, we also want such a solution to be rather generic which requires 
some sort of stability. This would be achieved if the solution contains only one unstable 
mode, and then the generic spacetime may first approach this solution and gradually leave 
it, recovering the present Friedmann universe, where we expect the higher order terms 
become irrelevant. We find that the solution N3 may give one possible candidate for such 
a model. We now discuss a new scenario obtained from this solution in the next section. 
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Table 2: The eigenvalues of linear perturbation equations for 7 = —1. We show two cases 
of /9: (3 — f3s and one for which we find an interesting scenario. 





solution 




h 


five eigenvalues {cFa \ a = 1, 2, • • • 5) 


-0.0052083 


N2 


0.45413 


0.45413 


( -4.5893, -4.5413, -3.6670, -0.87426, 0.048032 ) 


N3 


0.79802 


0.10781 


( -4.1012, -3.1488, -2.3823, -0.76653, 0.95268 ) 


N4 


0.50754 


0.43025 


( -4.5344, -4.4768, -3.6344, -0.89998, -0.057629 ) 


-0.2025 


Ni 


0.28195 


-0.39104 


(-10.2506, 0.9457 ± 1.8339Z, 1.8914, 12.1421) 


N2 


0.25005 


0.25005 


( -5.0572, -2.5005, -1.2502 ± 5.7723z, 2.5567 ) 


N3 


0.71567 


0.12395 


( -3.0268, -3.0147, -1.5073 ± 3.0912i, 0.01217 ) 


N4 


0.70803 


0.12661 


( -3.01038, -2.9986, -1.5052 ± 3.2449z, -0.011765 ) 



3 A Scenario for Large Extra Dimensions 

Let us discuss the evolution of the early universe for the solution N3(/i3, u^) = (0.79802, 0.10781) 
with j3 = j3s- The scale factor of the external space expands as e^-^^^*^^*. For a successful 
inflation (resolution of flatness and horizon problems), we need at least 60 e- foldings. Let 
us assume that inflation will end after 60 e-foldings, i.e. 0.79802tend ~ 60. The infla- 
tion will end because N3 has one unstable mode. During inflation the internal space also 
expands exponentially. When inflation ends, its scale becomes e'''^°''*^*<="'* pa 4000 times 
larger than the initial scale length, which we assume to be the IID Planck length (mj^/). 
After inflation, if the internal space settles down to static one, the present radius of extra 
dimensions is Rq ~ 4000mJ"/. Since this is slightly larger than the fundamental scale 
length, we may adopt the model of large extra dimensions, which was first proposed as a 
brane world by Arkani-Hamed et al. [15]. In this model, the 4D Planck mass is given by 

ml ~ Rlm\^ ~ 1.6 x lO^^m^^ . (14) 

We then find 

mil ~ 2.5 X 10"^^m4 ~ 600TeV . (15) 

This is our fundamental energy scale. The present scale of extra dimensions is 4000m]"j^^ ~ 
7 TeV~^, which could be observed in the accelerators of next generation. 
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We can also put our argument in a different way. Suppose that the e-folding of inflation 
is N, which is related to the stability of the solution. The 3-space expands as = e'^^"'"^ , 
while the internal space becomes e'^*'""^ times larger. It follows from Eq. (14) that 

mil ~ e^^^m4 . (16) 

Since mn > 1 TeV from the present experiments, we have a constraint on the e-folding as 
A*" < 10/x/ u. Then if we have TeV gravity and n > 6z/(> 0), we can naturally explain why 
the e-folding of inflation is about 60 and but not so large. Recall that the solution N3 

with 7 < gives 5.72 < fi/u < 10.22 (corresponding to 57< A^< 102) for any value of f3. 

(3) 

Although the above solution N3 has one unstable mode, its eigenvalue (T5 is of the 
same order of magnitude as other eigenvalues of stable modes as seen from Table 2 and 
is a little too large to give enough expansion. If the eigenvalue of the unstable mode is 
much smaller than those of other four stable modes, a preferable generalized de Sitter 
solution is naturally obtained for a wide range of initial conditions. Can we find such a 
possibility in superstrings or M-theory? We note that our starting Lagrangian has some 
ambiguity, that is, the fourth-order correction term 5*4 is fixed up to the Ricci curvature 
tensors. If we include correction terms including the Ricci curvature tensors, our basic 
equations will be modified. We might effectively take their effect into account by changing 
the value of our coefficient 7 or after the rescaling. Thus we may look for a preferable 
solution by changing /3. We find that the solution N3 with P = —0.2025 shows interesting 
behaviors (see Table 2). Four modes arc stable and the eigenvalue of one unstable mode 
is very small, i.e. = 0.01217. Then the time scale in which this unstable mode 

becomes important is evaluated as tus ~ (^'F'*)"^ ^ §2. Since the eigenvalues of other 
stable modes are of order unity, for a wide range of initial conditions, general solutions 
first approach the solution N3, which gives us an inflationary stage, after one dynamical 
time {t ~ 0(1)). The unstable mode becomes important at i ~ t^s, and then the inflation 
ends. We may have enough e-folding time of inflation {N /isfus ~ 58.8). In this case, 
we find mn ~ 3.3 x 10~-'^^m4 ~ 0.8 TeV, which gives us a TeV gravity theory. 

Although we find a successful exponential expansion and its natural end, this is not 
enough for a successful inflation. We need a reheating mechanism and have to create a 
density fluctuation as a seed of cosmic structure. A gravitational particle creation may 
provide a reheating mechanism [17], because the background spacetime is time dependent 
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and there might have some oscillation when the internal space settles down to static one 
which is required to explain our present universe. As for a density perturbation, our model 
may not give a good scenario because our energy scale is |7|~^''^ ~ Smu ~ 4 TeV. We 
have to invoke other mechanism for density perturbations such as a curvaton model [18] . 
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